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LONG-TERM  GOALS 

The  work  described  here  supports  the  Office  of  Naval  Research  Departmental  Research  Initiative 
(DRI)  for  research  on  the  “Impact  of  Typhoons  on  the  Western  Pacific  Ocean”  (ITOP).  The 
theme  of  the  DRI  is  to  better  characterize  and  predict  the  ocean  boundary  layer  (OBL)  and  its 
impact  on  typhoon  (hurricane)  evolution.  This  is  one  component  of  developing  improved  pre¬ 
diction  models  for  the  coupled  atmosphere-ocean-wave  system.  Cooling  of  the  sea  surface  tem¬ 
perature  (SST)  is  a  critical  coupling  variable  influencing  atmosphere-ocean  hurricane  dynamics; 
SST  is  largely  determined  by  OBL  turbulence,  surface  wave  processes,  and  mixed  layer  entrain¬ 
ment.  Our  research  goal  is  to  model  the  strongly  forced  wind  and  wave  driven  upper  OBL  using 
turbulence  resolving  large-eddy  simulation  (LES)  with  explicit  wave  effects,  viz.,  wave-current 
interactions  and  breaking  waves  and  examine  their  impact  on  ocean  mixing  during  hurricane 
events. 


OBJECTIVES 

The  specific  research  objectives  for  ITOP  are:  (1)  conduct  process  studies  using  LES  of  the 
OBL  with  different  combinations  of  time  varying  large  scale  forcings  and  surface  wave  effects 
and  examining  their  impact  on  ocean  mixing;  (2)  evaluate  and  compare  these  LES  results  with 
predictions  obtained  using  a  1-D  column  model  of  the  OBL  based  on  the  K-Prohle  Parameteri¬ 
zation  (KPP);  and  (3)  compare  our  simulation  results  with  available  observations.  Inertial  res¬ 
onance,  storm  residence  time,  and  the  larger  scale  environment  are  some  of  the  processes  to  be 
examined  in  our  simulations. 


APPROACH 

Present  computer  power  is  insufficient  to  simultaneously  resolve  all  the  dynamical  scales  in  a 
hurricane  driven  ocean  basin  where  the  largest  scales  of  motion  are  G(100s)  of  kilometers  in 
horizontal  directions  while  the  smallest  critical  scales  associated  with  spray  and  bubbles  in  the 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 
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air-sea  interface  are  0(1)  mm.  Our  interest  is  the  strongly  forced  wind  and  wave  driven  upper 
ocean  mixed  layer  where  the  turbulent  motions  are  0(1  —  500)  m.  In  order  to  examine  this 
range  of  scales  we  developed  an  LES  model  of  the  OBL  that  accounts  for  surface  wave  effects 
through  non-conservative  breaking  waves  and  phase-averaged,  conservative  wave-current  inter¬ 
actions  which  lead  to  Langmuir  circulations  (Sullivan  et  al.  2007).  The  equations  for  the  re¬ 
solved  flow  components  are  the  Craik-Leibovich  (CL)  theory  with  crucial  wave-current  coupling 
through  the  vortex  force.  The  larger  scale  forcing  (momentum  and  scalar  fluxes)  and  wave  fields 
are  externally  imposed  in  the  LES. 


WORK  COMPLETED 

Meetings:  We  attended  a  PI  meeting  in  Honolulu  HI  that  focused  on  refining  elements  of  the 
ITOP  held  campaign.  The  discussion  concentrated  on  the  deployment  of  the  numerous  drifters, 
buoys,  and  boats  and  the  coordination  with  the  aircraft.  A  summary  of  the  discussion  is  avail¬ 
able  from  E.  D’Asaro  (U.  Washington).  We  also  presented  a  summary  of  our  modeling  results 
from  the  HRES  and  ITOP  (Impacts  of  Typhoons  on  the  Western  Pacific  Ocean)  departmental 
research  initiatives  at  the  ONR  program  review  (June,  2009)  in  Chicago. 

Publication:  Last  year  we  wrote  an  article  for  Annual  Review  of  Fluid  Mechanics  (see  Sullivan 
&  McWilliams,  2010).  Our  contribution  discusses  the  coupling  processes  between  surface  gravity 
waves  and  adjacent  winds  and  currents  in  the  marine  boundary  layers.  Wind-wave  and  wave- 
current  interactions  are  important  as  these  processes  modulate  the  exchanges  of  momentum, 
heat,  and  gases  between  the  atmosphere  and  ocean.  Atmospheric  processes  we  discuss  are  wind- 
driven  waves,  wave-driven  winds,  steep  waves  and  drag  laws.  In  the  ocean  boundary  layer,  we 
summarize  the  currently  accepted  views  of  wave-current  interactions,  the  modeling  of  Langmuir 
circulations,  impacts  of  breaking  waves,  and  the  implications  of  wave-current  interactions  for 
mixing  at  high  winds. 

Modeling  Work:  During  this  Erst  phase  of  ITOP  we  concentrated  on  adapting  our  large- 
eddy  simulation  code  for  the  OBL  to  a  high  wind  regime  with  time  varying  winds,  surface  scalar 
fluxes,  and  wave  fields.  For  our  LES  process  studies,  Hurricane  Frances  is  selected  as  a  canoni¬ 
cal  storm.  Frances  was  a  large  category  4  hurricane  that  developed  in  the  Atlantic  basin  in  2004 
and  was  one  of  the  most  heavily  studied  storms  in  the  ONR-sponsored  Coupled  Boundary  Layer 
Air-Sea  Transfer  (CBLAST)  program  (Black  et  al.  2007).  Novel  profiling  boats  (Sanford  et  al. 
2007)  document  the  upper  ocean  mixing  induced  by  this  storm.  The  impact  of  the  storm  on 
currents  and  scalars  is  also  extensively  examined  by  Zedler  (2007)  using  the  MIT  ocean  model¬ 
ing  system.  Zedler  (2007)  estimates  Hurricane  Frances  had  a  radius  of  maximum  winds  R  =  40 
km  with  a  surface  rotational  wind  speed  Vr  ~  50ms-1. 

The  general  design  of  the  LES  experiments  is  to  pass  a  hurricane  vortex  over  small  LES  do¬ 
mains  as  sketched  in  hgure  1.  Then  the  surface  buxes  and  wave  fields  imposed  at  the  water  sur¬ 
face  vary  with  time  depending  on  the  storm  propagation  speed  Vs  and  the  location  of  the  LES 
domain  ( xleSiUles )  within  the  storm  track  (see  discussion  in  RESULTS).  Our  process  studies 
examine  the  impact  of  resonant  and  anti-resonant  wind  and  inertial-current  forcing  (e.g.,  Price, 
1981)  by  locating  LES  domains  to  the  right  and  left  of  the  storm  track.  Also,  the  solution  sen¬ 
sitivity  to  surface  cooling  is  considered,  i.e.,  simulations  with  and  without  a  surface  cooling  bux 
are  performed. 

The  initial  temperature  state  of  the  OBL  for  all  the  LES  experiments  is  neutrally  stratibed  over 
the  region  -33  m  <  z  <  0  m  bounded  by  a  stable  inversion  of  0.05  C  m-1  for  z  <  —33  m.  The 
LES  domain  (Xl,  Yf,  Zl)  =  (750,  750,  — 195)m  is  discretized  using  (Nx,  Ny,  Nz )  =  (500,  500, 160) 
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gridpoints;  the  vertical  spacing  varies  smoothly  from  Az  =  1.0  m  near  the  water  surface  to 
about  Az  =  1.4  m  near  the  lower  boundary.  Each  simulation  is  Erst  run  from  a  cold  start  for 
a  period  of  ~  10  hours  to  generate  turbulence.  The  LES  is  then  restarted  using  these  archived 
turbulent  flow  fields  as  initial  conditions.  For  LES  experiments  with  wave  effects,  we  set  the  tur¬ 
bulent  Langmuir  number  Lat  =  0.3  and  wave  age  =  0.6  which  Exes  the  parameters  in  an  ex¬ 
ponential  Stokes  drift  profile  (Sullivan  &  McWilliams,  2010).  As  shown  in  figure  2,  the  winds 
vary  considerably  with  t  and  as  a  result  the  computational  timestep  At  decreases  from  about 
7.1  seconds  at  the  beginning  of  the  simulation  to  about  0.65  seconds  at  the  time  of  maximum 
winds.  More  than  190,000  timesteps  are  needed  to  cover  the  entire  60  physical  hours  of  inter¬ 
est  for  cases  located  along  the  path  of  maximum  wind  speed.  We  note  that  the  OBL  LES  code 
used  here  is  a  variant  of  the  new  highly  parallel  atmospheric  code  that  is  being  developed  for 
the  high  resolution  air-sea  interaction  (HRES)  DRI.  The  algorithm  is  based  on  an  incompress¬ 
ible  Boussinesq  flow  model  and  solves  the  governing  equations  utilizing  a  mixed  finite-difference 
pseudospectral  scheme.  The  parallelization  is  accomplished  utilizing  the  Message  Passing  Inter¬ 
face  and  a  2-D  domain  decomposition,  see  Sullivan  &  Patton  (2008)  for  further  details. 


RESULTS 

Figure  2  shows  the  variable  time  forcing  applied  to  the  water  surface  on  the  right-hand-side  of 
the  storm  track  and  the  resulting  profiles  of  vertical  scalar  flux  and  vertical  velocity  variance. 

At  each  time,  the  profiles  of  ( w'6 )  and  (w12)  are  obtained  by  averaging  over  a  time  window  of 
about  30  minutes.  The  scalar  flux  and  vertical  velocity  variance  are  normalized  by  the  surface 
heat  flux  Q *  (shown  in  the  upper  panel  of  figure  2)  and  the  surface  friction  velocity  u2a  found 
from  a  saturated  bulk  aerodynamic  drag  rule  (Zedler,  2007),  respectively.  The  temporal  varia¬ 
tion  of  the  maximum  entrainment  flux  is  interesting.  Over  the  period  t  =  [0,  20]  hrs  the  max¬ 
imum  entrainment  flux  ratio  — ( w'e9)/Q *  ~  [0.2,  0.5],  i.e.,  the  entrainment  flux  is  only  slightly 
larger  than  is  observed  in  the  daytime  convective  boundary  layer  over  land  ( e.g Sullivan  et 
al.  1998).  We  attribute  this  increase  to  the  presence  of  Langmuir  cells  shown  in  figure  4  as  dis¬ 
cussed  in  McWilliams  et  al.  (1997).  For  t  >  20  hrs  the  rapid  increase  in  surface  winds  drives 
the  OBL  into  a  shear  dominated  regime  and  then  the  OBL  temperature  evolution  is  dominated 
by  entrainment  dynamics.  The  maximum  entrainment  flux  —  (w'e6) /Q*  ~  12  occurs  near  t  ~  35 
hrs,  i.e.,  at  the  time  of  maximum  winds.  Beyond  f  >  40  hrs  the  evolution  of  the  temperature 
field  is  complex  as  the  surface  wind  and  wave  forcing  relaxes  and  the  OBL  turbulence  decays  on 
a  background  of  strong  stable  stratification.  Also  the  hurricane  wake  becomes  strongly  three  di¬ 
mensional  (Sanford  et  al. 2007)  which  is  an  effect  not  captured  in  our  small  scale  LES  domains. 
The  vertical  velocity  variance  shows  a  layer  between  0  <  —  z  <  30  m  of  enhanced  vertical  veloc¬ 
ity  variance.  The  maximum  ( w'2)/u20  is  about  3.5  and  persists  over  the  majority  of  the  simula¬ 
tion  period  16  <  t  <  54  hrs.  This  elevated  vertical  velocity  variance  is  observed  in  all  LES  with 
vortex  force  under  steady  wind  forcing  with  the  amplitude  of  the  vertical  velocity  maximum  de¬ 
pendent  on  the  wave  field. 

The  impact  of  track  location  on  the  sea  surface  temperature  is  given  in  figure  3.  In  these  ide¬ 
alized  LES  the  cooling  is  greatest  on  the  resonant  right-hand-side  due  to  inertial  shear  insta¬ 
bilities,  as  expected.  The  inclusion  of  vortex  force,  based  on  a  simple  model  for  Stokes  drift, 
lowers  the  sea  surface  temperature  by  about  —(0.11,0.41,0.22)  degrees  C  or  (43,37, 15)%  at 
t  =  (30,40,50)  hours.  This  is  primarily  a  consequence  of  the  vigorous  Langmuir  circulations 
that  develop  in  the  mixed  layer  (see  figure  4).  Comparing  the  vertical  velocity  at  similar  depths 
we  notice  that  organized  streaky  structures  elongated  in  the  y— direction  form  in  the  presence  of 
vortex  force.  At  this  depth,  the  Langmuir  cells  have  length-to- width  aspect  ratios  of  about  8:1 
with  the  instantaneous  downwelling  speeds  in  these  structures  exceeding  -0.25  m  s-1. 
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IMPACT  /APPLICATIONS 


The  LES  results  obtained  here  for  hurricane  driven  OBLs  can  be  used  to  guide  the  interpreta¬ 
tion  of  observations  collected  during  the  ITOP  program.  In  addition,  the  results  can  be  used  to 
test  simpler  1-D  parameterizations  of  the  ocean  mixed  layer  that  are  used  in  large  scale  models. 
Our  particular  interest  is  in  evaluating  and  improving  the  so-called  K-profile  parameterization 
(KPP)  that  is  routinely  used  in  the  Regional  Ocean  Modeling  System  (ROMS). 


TRANSITIONS  &  RELATED  PROJECTS 

The  present  work  has  links  to  the  ONR  DRI  on  High  Resolution  Air-Sea  Interaction  (HRES) 
that  focuses  on  the  interaction  of  waves  and  turbulence  in  the  atmospheric  surface  layer.  The 
LES  model  being  developed  for  HRES  is  also  being  used  in  the  present  work. 
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Figure  1:  Sketch  of  an  idealized  hurricane  vortex  used  for  setting  the  surface  conditions  in  LES.  The 
vortex  is  propagating  upward  with  speed  Vs  and  has  a  characteristic  length  scale  (radius  of  maximum 
winds)  A.  The  families  of  solid  and  dotted  circles  indicate  the  position  of  the  vortex  at  initial  time  t 
and  at  a  later  time  t  +  At,  respectively.  The  fixed  LES  domains  to  the  right  (R)  and  left  (L)  of  the 
vortex  center  (c)  feel  the  time  history  of  wind  speed  and  direction  along  the  dotted  vertical  lines. 


6 


Figure  2:  Vertical  profiles  of  turbulent  scalar  flux  and  vertical  velocity  variance  from  an  LES 
of  an  OBL  driven  by  Hurricane  Frances  winds  and  buoyancy  with  vortex  force.  The  time  vari¬ 
ation  of  the  wind  (U,  V)  and  buoyancy  forcing  Q*  is  given  in  the  upper  panel.  The  LES  do¬ 
main  is  located  on  the  right  hand  (resonance)  side  of  the  storm,  initially  (55,  700)  km  east  and 
north  of  the  storm  center.  The  storm  propagates  northward  at  5.5  ms-1.  The  total  scalar  flux 
{w'9')/Q*  and  vertical  velocity  variance  ( w,2)/ul0 ,  are  shown  in  the  middle  and  lower  panels, 
respectively.  The  black  line  is  the  depth  of  the  maximum  vertical  temperature  gradient,  which 
can  be  interpreted  as  the  OBL  depth.  The  mesh  is  500  x  500  x  160  gridpoints  and  the  timestep 
At  varies  from  15  s  to  0.4  s  over  the  length  of  the  simulation.  The  total  number  of  timesteps 
~  200,000. 
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Figure  3:  Variation  of  SST  for  three  different  LES  driven  by  time  varying  winds  and  cooling, 
a)  resonant  side  of  storm  track  and  no  wave  effects;  b)  resonant  side  of  storm  track  with  simple 
model  of  Stokes  drift;  and  c)  non-resonant  side  of  storm  track  and  no  wave  effects. 
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Figure  4:  Visualization  of  the  vertical  velocity  field  in  an  x  —  y  plane  at  z  =  —30.0  m  at 

the  time  of  maximum  winds  shown  in  figure  2,  t  ~  35.3  hrs.  The  left  panel  shows  results  from 
an  LES  with  no  vortex  force  while  the  right  panel  is  a  simulation  with  Lat  =  0.3  and  wave  age 
=  0.6.  At  this  simulation  time,  the  wind  and  wave-propagation  directions  are  south  to  north.  In 
the  right  panel,  the  elongated  regions  with  w  «  0  (the  blue  contours)  are  the  downwelling  sig¬ 
natures  of  Langmuir  circulations.  The  color  bar  is  in  units  of  ms'1. 
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